/IDA071 4&0 


1 Ptf’7!-^  • 

I I,*  If  a*  i>T 


NONLINEAR  DYNAMIC  BUCKLING  OF 

THE 

SPHERICAL  CAPS  WITH 

' 

GEORGE 

INITIAL  IMPERFECTIONS 

P 

X 

WASHINGTON 

UNIVERSITY 

-f 

f >- 

Q_ 

O 

C_? 


1 


STUDENTS  FACULTY  STUDY  R 
ESEARCH  DEVELOPMENT  FUf 
URE  CAREER  CREATIVITY  CC 
MMUNITY  LEADERS 
NOLOGY  FRONT 
ENGINEERING  AP 
GEORGE  WASHI 


LIU' 


Lie.. 


SCHOOL  or  LNC.INEERING 
AND  APPLIED  SCIENCE 


THIS  I 


: RELEASE  AND  SALE:  ITS  DISTRIBUTION  IS  UNLIMITED. 


NONLINEAR  DYNAMIC  BUCKLING  OF 
SPHERICAL  CAPS  WITH 
INITIAL  IMPERFECTIONS 

I ^ 


Robert  Kao 


Sponsored  by 

Office  of  Naval  Research 
Arlington,  Virginia  22217 


Contract  Number 


NAVY  00014-75-C-0946 


June  1979 


School  of  Engineering  and  Applied  Science 
The  George  Washington  University 
Washington,  D.C.  20052 


Jt) 


'%  9 0 a 


i 'O 


\J  ** 


NONLINEAR  DYNAMIC  BUCKLING  OF 
SPHERICAL  CAPS  WITH  INITIAL  IMPERFECTIONS1 


By 

Robert  Kao 


The  research  reported  on  here  was  supported  by  the  Office 
of  Naval  Research,  Contract  Number  NAVY  00014- 75-C- 0946 


ABSTRACT 


A finite  difference  method  is  developed  for  the  large  de- 
formation elastic-plastic  dynamic  buckling  analysis  of  axisym- 
metric  spherical  caps  with  initial  imperfections.  The  problem 
formulation  is  based  on  governing  differential  equations  of 
motion,  treating  the  plastic  deformation  as  an  effective  plas- 
tic load.  Both  perfectly  plastic  and  strain  hardening  behavior 
are  implemented  in  the  program.  Strain  hardening  is  incorporated 
through  use  of  the  Prager- Z iegler  kinematic  hardening  rule,  so 
that  the  Bauschinger  effect  is  accounted  for.  The  solution  for 
the  large  deformation  elastic-plastic  dynamic  response  of  a 
spherical  cap  is  compared  very  favorably  with  other  findings. 

Two  spherical  cap  models  are  selected  to  study  the  title  problem. 
Results  obtained  indicate  that  both  plastic  yielding  and  initial 
imperfection  play  significant  roles  in  reducing  the  load  carry- 
ing capacity  of  these  shell  structures.  Both  increase  their 
influence  as  the  thickness  to  radius  ratio  and  the  imperfection 
magnitude  increase,  respectively.  It  is  also  found  that  dynamic 
effect  has  the  influence  of  lowering  load  carrying  capacity  of 
perfect  spherical  caps;  however,  its  influence  on  imperfect 
spherical  caps  depends  on  the  magnitude  of  initial  imperfections. 


INTRODUCTION 


Dynamic  buckling  analysis  of  shell  structures  has  been  the 
subject  of  intense  study.  Shell  structures  designed  according 
to  quasi-static  analysis  may  fail  under  conditions  of  dynamic 
loading.  For  a more  realistic  prediction  on  the  load  carrying 
capacity  of  these  shells,  in  addition  to  the  dynamic  effect, 
considerations  should  also  include  other  factors  such  as  non- 
linearities  in  both  material  and  geometry,  initial  imperfec- 
tions, etc.,  since  these  factors,  in  a different  manner,  may 
also  affect  the  magnitude  of  this  capacity. 

Large  deformation  dynamic  buckling  analysis  was  studied 
by  several  authors,  notably  Budiansky  and  Roth  [1],  Huang  [2], 
Stephens  and  Fulton  [3],  Stricklin,  et  al.  [4]  and  Ball  and 
Burt  [5].  A striking  difference  between  the  static  and  dynamic 
buckling  analyses  is  on  the  buckling  criterion.  The  dynamic 
snap-through , suggested  in  [1]  and  adopted  in  all  these  studies, 
is  generally  accepted  as  a means  to  obtain  dynamic  buckling 
loads  of  spherical  caps  under  uniform  step  loading.  Results  of 
spherical  caps  obtained  from  these  studies  are  in  reasonable 
agreement . 

Lffect  of  initial  imperfections  on  dynamic  buckling  loads 
of  axisymmetric  spherical  caps  was  examined  in  Ref.  [6].  It  is 
found  that  initial  imperfections  do  indeed  have  the  effect  of 
reducing  the  cap  buckling  capacity,  and  that  the  rate  of  change 
of  the  buckling  load  with  respect  to  the  imperfection  is 
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greater  when  caps  respond  statically  than  dynamically.  Two 
types  of  loading  were  also  treated  in  Ref.  [6],  namely,  step 
pressure  of  infinite  duration  and  triangular  pulses  of  various 
time  durations.  It  is  also  revealed  that  the  pulse  duration 
lias  a very  profound  impact  on  the  magnitude  of  the  dynamic 
buckling  load,  and  that  the  step  loading  of  infinite  duration, 
as  a limiting  case  of  triangular  pulses,  provides  the  most  severe 
loading  situation  for  the  dynamic  buckling  analysis. 

The  combined  effect  of  material  and  geometric  nonlineari- 
tics  on  the  shell  behavior  was  studied  by  Marcal  [7],  Bushnell 
[81,  Levine,  et  al.  [9],  and  Kao  [10];  the  effect  of  initial 
imperfections  was  also  considered  in  Ref.  [10].  It  is  observed 
in  Ref.  [10]  that  both  initial  imperfection  and  plastic  de- 
formation have  a similar  influence  of  reducing  shell  load 
carrying  capacity,  and  that  the  influence  of  plastic  de- 
formation decreases  with  increase  of  the  imperfection  magnitude. 
It  is  also  found  from  this  reference  that  the  cap  geometric 
parameter  \ , an  important  factor  in  the  elastic  response,  be- 
comes meaningless  in  the  elastic-plastic  buckling  analysis  of 
spherical  caps. 

The  purpose  of  this  paper  is  to  obtain  large  deformation 
clastic-plastic  dynamic  buckling  loads  of  axisymmetric  spherical 
caps  with  initial  imperfections.  The  problem  formulation  is 
based  on  equilibrium  equations  of  motion,  treating  the  plastic 
deformation  as  an  effective  plastic  load  [10,11].  Governing 
equations  are  replaced  by  finite  differences  and  a direct 
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method  Ls  applied  for  the  time  integration.  At  each  time  step, 
a relaxation  technique  [12]  is  introduced  to  solve  nonlinear 
algebraic  equations.  Plasticity  relations  are  derived  from  the 
incremental  flow  theory  [13]  together  with  the  von  Mises  in- 
itial yield  criterion  and  the  Prager- Z iegler  kinematic  harden- 
ing rule  [14,15];  this  rule  predicts  an  ideal  Bauschinger 
effect . 

In  the  next  two  sections,  governing  equations  of  motion 
and  plasticity  relations  are  subsequently  given.  This  is  fol- 
lowed by  the  description  of  the  general  solution  procedure  and 
dynamic  buckling  criteria.  In  the  NUMERICAL  RESULTS  AND  DIS- 
CUSSION' section,  first,  the  present  solution  of  the  large  de- 
formation elastic-plastic  dynamic  response  of  a spherical  cap 
is  compared  with  other  findings;  then,  results  of  the  title 
problem  are  presented  and  a genera]  discussion  on  these  re- 
sults is  given.  Conclusions  are  outlined  in  the  final  section. 

GOVERN  1 NG  EQUATIONS  OF  MOTION 

The  geometry  of  a spherical  cap  is  shown  in  Fig.  la,  in 
which  II  is  the  central  height,  R the  shell  radius  and  a is  the 
base  radius;  W(r,t)  and  U(r,t)  are  the  displacements  along  nor- 
mal and  tangential  directions,  respectively,  and  hA(r,0)  is 
the  initial  imperfection;  q(r,t)  is  the  applied  pressure.  The 
undeformed  shape  of  the  perfect  spherical  cap  can  be  adequately 
described  by 


Z = H [1  - (r /a)2] 


(la) 
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where  r is  the  radial  coordinate.  The  radius  of  curvature  of 
the  shell  is  approximated  by 

R = a2/  211  (lb) 

Figure  lb  shows  an  element  with  the  stress  resultants  Nr 

and  N„,  the  transverse  shear  Q and  the  moments  M and  ML . The 
0 *r  r 9 

equations  of  equilibrium  including  inertia  forces  in  this  ele- 
ment are 


(rNr)  - Nq  = 0 

[rNr  (Wf  - Z)  + rQr]  + rq  = rphW 
(rMr)  - M0  - rQr  = 0 


(2) 

(3) 

(4) 


where  h is  the  shell  thickness,  p is  the  mass  per  unit  volume 
of  the  shell  and  = IV  + VL  ; prime  and  dot  denote  differen- 
tiations with  respect  to  r and  t,  respectively,  and  0 is  the 
circumferential  coordinate. 

It  is  noted  that,  in  light  of  the  assumed  shallowness, 
the  effect  of  transverse  shear  and  radial  inertial  force 
rphU  in  Eq.  (2)  is  neglected,  and  that  the  nonlinearity  has 
been  introduced  in  Eq.  (3)  by  considering  the  influence  of  W^. 

Eliminating  Qr  in  Eqs.  (3)  and  (4),  we  obtain 


M'  ' 
r 


2 

- M' 
r r 


- - Ml 
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Nr(W” 


Wf 

N.  ( — 
0 


V* 


+ q = phW  .(5) 


Eqs.  (2)  and  (5)  are  the  basic  equations  for  the  dynamic 
analysis  of  axisymmetric  spherical  caps  in  this  paper.  To 
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solve  these  two  equations,  we  elect  to  express  them  in  terms 
of  the  components  U and  W of  the  displacement.  For  this  pur- 
pose, stress- strain  and  strain-displacement  relations  are  re- 
quired. 


Stress-Strain  and  Strain-Displacment  Relations 


For  a shell  deforming  into  a plastic  range,  the  strain  in 
a point  within  the  thickness  can  be  considered  as  a combination 
of  its  elastic  and  plastic  components: 

(e)  = { ee } + {ep}  (6) 

where  {e},  { e° } and  {ep}  are  the  total,  elastic  and  plastic 
strain  vectors,  respectively.  The  total  strain  can  also  be  ex- 
pressed as  a sum  of  the  membrane  and  bending  components: 

{el  = { e } + z{k}  (7) 

where  z is  the  vertical  coordinate  through  the  shell  thickness 
(Fig.  lb). 

Membrane  and  bending  strains  are  related  to  displacements 
by 


Er  = U'  - | + 1(W’)2  + W'W! 

U W 
L 0 r ' R 

Kr  = -W" 
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Note  that  the  elastic  components  of  strains  are  the  only 
strains  which  can  be  related  to  stresses  by  Hook's  law: 

{a}  = [E]  ({e}  - {ep})  (9) 

where  [E] , the  elastic  strain  to  stress  transformation  matrix, 
is  given  as 


[E]  = 


M' ;] 


(10) 


in  which  E is  Young's  modulus  and  v is  Poisson's  ratio. 

Membrane  stress  resultants  and  bending  moments  are  ob- 
tained by 


h/2 

(N)  = I (o)dz 
-h/2 


(11) 


(M)  = 


/ 


h/2 


{ a } zdz 


(12) 


-h/2 


Substituting  Eqs.  (6-9)  into  Eqs.  (11)  and  (12),  we  obtain 
the  membrane  forces 


Eh  P 

1-v2  Lv  1. 


where  the  effective  plastic  membrane  forces  are 


dz 


Nj 

r h/2 

jep  ) 

' n 

= fE]  / 

r i 

N?  I 

" h / 2 

ee  \ 

(13) 


(14) 
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and  the  moments 
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where  D = Eh3/12(l-v2)  , 


and  effective  plastic  moments  are 


zdz 


Governing  Equations  in  Displacements 


(15) 


(16) 


The  governing  equation  involving  the  first  major  displace- 
ment U can  now  be  obtained  by  substituting  Eq.  (13)  into  Eq.  (^) : 

U"  + T~  ' ^2  + G(W)  = TF  ql  (17) 

where 

G (W)  = F^(W)  + vF’(W)  + (Fr-F0) 

Fr(W)  = -W/R  + ( W * ) 2 / 2 + W'W! 

F0(W)  = -W/R  (18) 

F'fW)  = -W'/R  + W'W"  + W’W'.'  + W'W! 

r v ' i i 

F ' (W)  = -W'/R 

and  c|‘j  , an  effective  plastic  load,  is  related  to  effective 
plastic  stress  resultants  by 


q 


P = 
1 


(nJ)’  + NlP/r-NP/r 


(19) 


y 


i 
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The  second  equation  involving  IV  as  its  major  displacement 
can  also  be  obtained  by  substituting  Eq.  (15)  into  Eq.  (5): 

Dv4W  - — — — 2-(e  +V£g)  (W'i+l/R)  - -^(e.+vel  (W> /r  + 1/R) 

1-v  r 0 1 l-vZ  0 r f 

= q - qP  - qP  - phW  (20) 

where  V4  = V^(V^)  and  ()=()"+  ( )'/r;  the  membrane 
strains  er,CQ  are  defined  in  Eqs.  (8),  effective  plastic  loads 
qP  and  qP  are  given  as 


qP  = NP(W”+1/R)  + NP  (W£/ r+l/R) 

(21) 

qP  = (MP)"  + 2 (MP) ' / r - (Mq) ’ /r 

(22) 

Equations  (17)  and  (20)  are  two  fundamental  governing 
equations  in  terms  of  displacements  for  the  present  analysis. 

Boundary  Conditions 

At  the  shell  apex,  the  nature  of  axisymmetry  requires  that 

W (0)  = 0 (23) 

U(0)  = 0 (24) 

Along  the  outer  edge  (r  = a) , if  the  cap  is  clamped: 

U(a)  = W(a)  = W (a)  = 0 . (25) 


On  the  other  hand,  if  the  cap  is  simply  supported,  it  re- 
quires that 
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U(a)  = W (a)  = 0 (26a) 

and  tliat  Mr(a)  in  Eq.  (15)  to  be  zero,  i.e. 

D(~I  + 7 37)  = -Mj(a)  , r = a (26b) 

where  is  defined  in  Eq.  (16). 


Nondimensional  Form 


For  convenience,  following  nondimensional  quantities  are 
int roduced : 


x = r /a 

7 7 7 

A = nr  a /Rh 

( )'  = 3(  )/3x 
(’)  = 3(  )/3t 

t = /l:/pR2  t 
wi  = IV. /h 


m4  = 12(l-v2) 
qcr  = 4Fh2/R2m2 

p(r,t)  = q(r,t)/qcr  (27) 
u = aU/h2 

w = W / h 


where  q^r  is  the  classical  buckling  pressure  of  a complete 
spherical  shell  of  the  same  radius  of  curvature  and  thickness. 

By  using  Fq.  (27),  nondimensional  forms  of  Fqs.  (17)  and 
(20)  become 

u"  + -2  + R(w)  = q?  (28) 

X x4  Eh3  1 


and 
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V"  w - 12(er  + ve  ) ( W’J-+  ^ 1 - 12(cQ  + ver)  I ^ \ 

m / \ m 


X 


4" 


» /\  IT1  <3  /-  P , P s -x  *+ 

= 4 -f  P T~(qZ qP  ' A w 


4 4 
1 c 

m~  El/ 

where  g(w)  and  f(w)  terms  are  nondimens ional  counterparts  of 
similar  terms  in  Eq.  (18) : 


(29) 


g(w)  = fj.(w)  + vf^(w)  + (1-v)  [fr  (w) 


f0(w)]/x 


fr(w) 


m 


L 17 

'~2  w + y(w')  + w'w!^ 


f9(w)  2 

m 


w 


(30) 


f;(w) 


m 


^5-  w'  + w'w"  + w'w'.'  + w"w! 
2 li 


f0(w)  " ' "T  w' 
m 


and  e , c„  are  nond imens ional  quantities  of  membrane  strains 
r 0 

e and  in  Eq.  (8) 


cr  = u'  - w + j(w')^  + w'wV 
in 


— u X 

eG  ' x ' 2 W 

m 


(31) 


CONSTITUTIVE  EQUATIONS  OF  PLASTICITY 


The  response  of  an  elastic-plastic  material  can  be  described 


[151  by 
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(a)  an  initial  yield  condition,  specifying  the  state 
of  stress  for  which  plastic  flow  first  sets  in, 

(b)  a flow  rule,  relating  the  plastic  strain  increment 
with  the  stress  and  the  stress  increment,  and 

(c)  a hardening  rule,  specifying  the  subsequent  yield 
condition  in  the  course  of  plastic  flow. 

In  this  paper,  the  von  Mises  yield  condition,  which  describes  a 
smooth  surface  in  stress  space  and  represents  a simple  mathema- 
tical function,  is  chosen  as  the  initial  yield  condition.  The 
flow  theory  of  von  Mises  and  the  Prager- Ziegler  kinematic  har- 
dening rule  are  also  selected;  this  rule  predicts  an  ideal 
Bauschinger  effect. 

In  a 9-space  stress  field  with  origin  0 (Fig.  2),  the  von 
Mises  initial  yield  surface  can  be  described  by 

2 

F(Oij)  .=  k = constant  (32) 

For  an  initially  isotropic  material,  the  form  of  the  function  F 
is  invariant  with  respect  to  a rotation  of  the  stress  state. 

The  Prager' s hardening  rule  [14]  assumes  that  during  plastic 
deformation,  the  yield  surface  moves  in  a translation  without 
changing  its  shape  and  hence  the  subsequent  yi^ld  surface  takes 
the  form 

F(°ij  - «ij)  = k2  (33) 

where  represents  the  total  translation  of  the  yield  surface 
center  which  is  a measure  of  the  degree  of  work  hardening.  In 
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the  space  Ojj>  °^j  is  the  position  vector  of  the  yield  surface 
center  C which  before  plastic  deformation  takes  place  is  located 
at  the  origin  (Fig.  2). 

The  flow  rule  of  von  Mises  gives  the  following  expression: 


deij  = 97T7  dA  , dX  > 0 (34) 

which  indicates  that  the  plastic  strain  increment  de? . lies  in 

1.1 

the  exterior  normal  of  the  yield  surface  (33) . 

The  Ziegler’s  modification  [15]  of  Prager's  hardening 
rule  suggests  that  the  surface  (33)  moves  in  the  direction  of 
the  radius  connecting  its  center  with  the  stress  point  (Fig.  2)  : 


daij  = (°ij "“ij )dF  > dp  > 0 (35) 

d X in  Eq.  (34)  and  dp  in  Eq.  (35)  are  to  be  determined. 

Determination  of  dp  is  based  on  the  condition  that  stress 
point  always  remains  on  the  yield  surface  in  plastic  flow. 

This  condition,  in  fact,  states  that  for  an  infinitesimal  in- 
crement of  loading  the  vector  do^  - dct^  must  be  orthogonal  to 
the  outer  normal  to  the  yield  surface: 


-ToT. 


3f(°irZAi>  _ 


Substituting  (35)  into  (36)  yields 


(36) 


dp 


(omn7%n]3F/3omn 


(37a) 


and  hence 
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aij  (a  -a  ) 3 F/ 30  ^°ii  aip 
J v ran  mrr  ' mn  J J 


(37b) 


According  to  Refs.  [14,15],  9A  in  (34)  can  be  obtained 
by  assuming  doi^  = cde'^  and  putting  this  relation  into  (36)  : 

„ i (3F/3t,ki)doki 

- c OF/So^HSF/Sc^) 


(38a) 


and  hence 


j P 1 (3F/3aki.)dokJt 

ij  " c (3F/9amn)(9F/ 


9F/9a  ) 9a.  . 

mn-*  ij 


(38b) 


where  c is  the  hardening  coefficient  and  can  be  determined 
from  the  uniaxial  stress-strain  relation. 

Now,  let  us  write  Eq.  (33)  in  terms  of  3 principal 
stresses : 


f = F - k‘ 


= \[(a\-a2^2  + (a2'°3)2  + ^a3'°l)2l  ' °y  = 0 

where  is  the  yield  stress  in  uniaxial  tension  and  a^  = 

°1  ' al»  a2  = °2  a2 ’ °3  = °3  ' a3  ’ 

For  the  case  of  plane  stress  (Fig.  3) , we  have 

°3  = a3  = °3  = 0 
and  Eq.  (39)  reduces  to 

—2  — — —2  0 
f = °1  ' °l  °2  + °2  ‘ = 0 

By  putting  Eq.  (39)  into  Eq.  (38b)  and  using  relations 
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(40-41) , we  obtain 


where 

_ (al~CT2^2^°y  > ^2  = ^ °2 " ° 1 1 ^ ) /°y  ’ ^ = 7 c 
From  Hook's  law,  we  have 


Introducing  (42)  into  (44)  provides 


where 


C11  * I <D*E  S2> 


c i 2 ' n CDv-nSls2)  - c21 


C22  = I <D*E  Sl> 

« = D(l-v2)  ♦ E(S2  + 2vSrS2  + S“) 


(42) 


(43) 


(44) 


(45) 


(46) 
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For  a given  displacement  field,  (Ae)  are  obtained  from  incre- 
mental forms  of  Eqs.  (7)  and  (3).  Eqs.  (46)  are  the  stress- 
strain  relations  during  the  course  of  plastic  flow  (loading). 
Otherwise,  Eqs.  (44)  should  be  used  for  all  other  stress  com- 
putat ions . 

Up  to  this  point,  the  hardening  coefficient  c in  Eq.  (43) 
is  the  only  constant  yet  to  be  decided.  If  the  structure  is  in 
a state  of  uniaxial  stress,  the  stress -strain  relation  is  the 
same  as  that  obtained  from  tension  or  compression  tests.  The 
expression  for  the  hardening  coefficient  c can  be  readily  ob- 
tained from  Eq.  (42)  by  setting  Aep  = Aep,  Aa  = Aa^: 


I = 3 Aep 

c 1 Aa 


(47a) 


or 

D = . (47b) 

Ac1 

This  equation  shows  that  the  value  of  D is  equal  to  the  slope 
of  the  uniaxial  stress-plastic  strain  curve. 

From  Eqs.  (47),  the  values  of  the  hardening  coefficient 
for  two  special  cases  may  be  specified  here:  (i)  for  an 
clast  ic- ideally  plastic  material,  c = D = 0,  (ii)  for  a linear 
hardening  material,  D = El;^/ (E-F.^ ) , where  E^  is  the  tangent 
modulus.  Uniaxial  st rcss- st ra in  curves  for  these  two  types 
of  hardening  material  are  given  in  Fig.  4. 
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Nonlinear  Hardening 


Also  given  in  Fig.  4 is  the  stress- strain  curve  for  a 
nonlinear  hardening  material.  Its  D value,  instead  of  being 
a constant,  depends  on  the  state  of  stress.  One  way  of  deal- 
ing with  this  rather  complicate  situation  is  based  on  the 
Rambcrg-Osgood  representation  of  a uniaxial  stress-strain 
curve  [ 18] : 


e = -rv  + 


where  n - 1 ♦ }°8-P7/7), . , 

log  (Vt'Vss* 

e is  the  total  strain,  E is  the  slope  of  the  linear  portion 
of  the  stress  - strain  curve,  and  ^ and  are  the  stress- 

es at  which  the  curve  has  secant  moduli  of  0.7E  and  0.85E, 
respectively. 

It  is  understood  that  the  nonlinear  term  in  Eq.  (48)  is 
the  plastic  strain.  D value  of  this  material  can  be  obtained 


(49a) 


We  may  generalize  this  equation  to  a multiaxial  state  of 
stress  [16,17]: 


D - 7E 

3n  — 

a 

where  a,  the  effective  stress,  is  defined  as 


(49b) 
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0 ■ °ij  °ij  (50a) 

Ln  which  ajj  is  the  deviatoric  stress.  For  the  case  of  plane 
stress, 

- 2 2 

° = °1  " °1  °2  + °2  * (50b) 

Loading  Criteria 

In  the  incremental  solution  procedure  for  elastic-plastic 
problems,  in  addition  to  the  constitutive  relations,  it  is 
necessary  to  have  a loading/unloading  criterion.  For  this 
purpose,  let's  define  f = (3f/9ai ^ ) dai ^ , where  f has  been 
defined  in  Fq.  (39).  Loading,  unloading  and  neutral  loading 

are  associated  with  the  plastic  state  f = 0,  and  are 

• • • 

characterized  by  f > 0 , f < 0 and  f = 0,  respectively.  When 
loading  or  neutral  loading  takes  place,  Eq.  (46)  must  be 
applied.  On  the  other  hand,  Eq.  (44)  must  be  used  when  un- 
loading occurs. 

■SOLUTION  STRATEGY 

For  convenience,  a simple  flow  chart  is  sketched  in  Fig.  5 
to  demonstrate  the  general  solution  procedure.  The  entire  pro- 
cess, which  is  to  obtain  for  a given  applied  load  q (or  p)  an 
elastic-plastic  transient  response  of  axisymmetric  spherical 
caps,  is  divided  into  two  major  loops,  namely,  the  elastic 
solution  and  material  property  loops. 
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In  the  elastic  solution  loop,  all  material  properties  at 
a specific  time  x are  held  constant.  Consequently,  the  effec- 
tive plastic  loads,  q'j*  in  governing  equations  (28)  and  (29) 
are  fixed  and  combined  with  the  actual  externally  applied  load 
q,  the  problem  is  thus  reduced  to  an  elastic  large  deformation 
prob 1cm . 

The  second  time  derivative  of  w in  F.q.  (29)  is  approximated 
by  the  Houbolt's  third-order  backwards  difference  expression  [19] 

w (x,x)  = ( 1 / 6 ^ ) [2w  ( x , x ) - 5w  (x  , x - 6) 

+ 4w  (x,x  - 26)  - w (x,x  - 36)]  (51) 

where  6 = Ax  is  the  equal  time  increment.  The  accuracy  of  this 

2 

representation  is  of  order  6 . Special  attention  is  devoted  to 
the  first  few  time  steps  where  Eq.  (51)  cannot  be  applied 
directly.  Before  giving  the  expression  for  these  first  few 
time  steps,  we  note  that  the  initial  conditions  are  of  the  form 

w (x , 0)  = 0,  w (x , 0)  = 0 (52) 

Prom  Eq.  (52),  expression  of  Eq.  (51)  for  the  first  few 
time  steps  can  readily  be  obtained  [2]  as 

1)  x = 0,  w (x,0)  = 0 

2)  x = 6;  since  w (x,0)  = 0,  we  have  w (x,-6)  = w (x,6) 

and  hence  w (x,6)  = (2/62)  w (x,6) 

3)  x = 26;  w (w , 26 ) = (2/62)  [w  (x,26)  - 3w  (x,6)] 

4)  x > 26,  Eq.  (9)  can  be  applied  directly. 
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At  this  point,  central  finite  differences  are  used  to  trans- 
form Eqs.  (28)  and  (29)  to  a discrete  system  of  equations,  and 
the  nonlinear  relaxation  technique  [12]  is  employed  to  solve 
these  nonlinear  equations  The  iteration  in  this  loop  is  con- 
sidered converged  when  the  average  absolute  change  of  u and  w 

displacements  at  all  points  is  less  than  .0001.  / 

k f 

With  the  new  displacement  field  (u)  , the  material  properties 

must  be  updated  so  that  the  (nonlinear)  stress  - strain  relation 
can  be  satisfied  at  all  points  over  the  shell  surface  and  through 
the  thickness.  In  the  material  property  loop,  the  loading 
criterion  is  first  checked.  If  the  material  point  is  in  an  un- 
loading situation  or  still  in  the  elastic  range,  Hook's  law  is 
used  and  the  computation  is  very  straight-forward. 

If  the  material  point  is  in  a loading  situation,  incremental 

strains  (Ae)  are  computed  from  F.q.  (7)  by  using  (Au)  = 
k k - 1 

(u)  - (u)  (where  k is  the  number  of  the  material  property 

loop,  see  Fig.  5),  and  (Ao)  are  computed  from  Eq . (46).  Having 
obtained  (Ao),  we  find  new  values  of  (o),  (ee)  and  {e^}.  The 
effective  plastic  loads  q?  are  then  evaluated  from  integration 
formulae  (19),  (21)  and  (22)  through  the  use  of  a Simpson's 

rule  (9  thickness  points  are  used  for  the  entire  computation  , 

4 

| 

in  this  paper) . 

The  material  property  is  considered  to  be  updated  if  aver- 
age absolute  change  of  displacements  between  the  present  and 
previous  material  property  loop  is  less  than  0.0005.  Otherwise, 
the  iteration  goes  to  the  elastic  solution  loop  and  the  entire 


operation  is  repeated  until  the  material  property  and  equilibrium 
equations  are  simultaneously  satisfied. 

In  our  numerical  computation,  the  number  of  nodal  points 
are  selected  such  that  a subsequent  increase  in  nodal  points 
does  not  significantly  affect  the  magnitude  of  the  static 
buckling  load.  With  this  consideration,  14  nodal  points  for 
X = 5 and  up  to  22  points  for  X = 10  are  adopted. 

A time  step  must  also  be  selected  very  judiciously.  A good 
selection  is  such  that  the  results  are  within  a desired  accuracy, 
but  not  too  small  in  light  of  computer  time  considerations.  A 
reasonable  compromise  of  equal  time  increments  of  6 = 0.10  is 
used;  this  same  time  increment  is  also  selected  in  Ref.  [3], 

The  axisyminetric  initial  imperfection  adopted  in  this 
study  is  of  the  dimple  type  which  was  also  used  in  Ref.  [6]. 

This  type  of  imperfection  is  expressed  mathematically  as 

wi  = (Wi0/h)d  ' x2)3  (53) 

where  W.()  is  the  maximum  imperfection  which  occurs  at  the  shell 
apex.  Selection  of  this  expression  is,  in  fact,  quite  arbi- 
trary. However,  it  does  provide  an  adequate  description  for 
actual  shells  since  the  important  parameter  is  the  maximum 
eccentricity  and  not  the  imperfection  shape  function. 

DYNAMIC  BUCKLING  CRITERION 

Criteria  for  dynamic  axisymmetric  buckling  of  the  shallow 
spherical  shell  arc  not  as  well  defined  as  for  static  buckling, 
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and  require  an  evaluation  of  the  transient  response  of  the  shell 
for  various  load  levels. 

The  criterion  adopted  most  widely  (Refs.  [1-6])  is  based 
on  plots  of  the  peak  nond imens iona 1 average  displacement  in  time 
history,  A , of  the  shell  structure  versus  the  amplitude  of 
the  load  where  A is  the  average  displacement  and  has  been  defined 
[3,5,6]  in  dimensionless  form  as  follows: 


A = 


G 

/ 


rWdr  / 


rZdr 


(54) 


The  numerator  is  the  volume  generated  by  the  shell  deformation 
and  the  denominator  is  the  constant  volume  under  the  cap. 

For  loads  below  the  dynamic  buckling  load,  where  the  non- 
linearities  are  small,  the  relationship  between  A and  the 
load  amplitude  is  approximately  linear.  However,  there  may 
exist  a certain  value  of  load  amplitude  at  which  a very  small 
increase  in  this  amplitude  produces  a very  large  in  A . This 
load  is  defined  as  the  critical  pressure  for  axisymmetric 
buckl ing. 


NUMFRI CAL  RESULTS  AND  DISCUSSION 

In  reference  [10],  a computer  program  was  successfully 
developed  for  the  large  deformation  elastic-plastic  analysis 
of  spherical  caps  with  initial  imperfections.  This  program  is 
extended  to  include  the  dynamic  effect  for  the  calculation  in 


22 


this  paper.  Before  proceeding  to  obtain  numerical  results  for 
the  title  problem,  it  is  advisable  to  solve  an  example  problem 
and  compare  the  present  result  with  other  existing  solutions. 

The  comparison  should  be  very  meaningful  if  other  solutions 
arc  obtained  using  complete  different  method  and  theory  from 
the  approach  adopted  herein. 

Comparison  with  Other  Solutions 

The  example  problem  considered  here  is  the  dynamic  response 

of  a shallow  spherical  shell  subjected  to  a uniformly  distributed 

2 

external  step  pressure  q = 600  lb/in.  Geometric  dimensions 
and  material  properties  are  outlined  in  Fig.  6.  A total  of  14 
stations  along  radial  direction,  9 thickness  points  and  a time 
step  of  At  = 10  ^ seconds  are  used  for  the  calculation.  The 
material  is  assumed  to  obey  von  Mises  initial  yield  criterion 
with  linear  kinematic  hardening. 

The  large  deformation  elastic-plastic  dynamic  response  is 
presented  in  Fig.  6,  which  also  shows  the  result  obtained  by 
assuming  the  cap  to  remain  clastic  throughout  the  analysis.  The 
static,  elastic-plastic  solution  is  also  given  in  this  figure. 

The  effect  of  material  nonlinearity  on  the  dynamic  behavior 
of  the  cap  is  significant.  It  is  observed  that  the  mean  value 
about  which  the  apex  displacement  oscillates  is  greater  in  the 
clastic-plastic  situation  compared  to  the  elastic  analysis 


wherein  the  vibration  occurs  around  the  static  displacement. 

Figure  6 also  displays  the  solutions  obtained  by  Nagarajan 
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I 

and  Popov  [20|  and  Bathe,  R a it> m and  Wilson  [21].  Both  solutions 
were  based  on  the  finite  element  formulation,  von  Mises  yield 
condition  with  linear  isotropic  hardening  and  a Newmark  gen- 
er  ized  acceleration  scheme.  Despite  differences  in  the  solu- 
tion method  and  theory  employed,  the  comparison  among  these 
results  is  excellent.  It  is  noted  that  the  same  time  step 
(10  ‘’sec.)  is  also  used  in  their  computations. 

With  the  validity  of  the  computer  program  established  for 
nonlinear  dynamic  analysis  as  demonstrated  in  this  example,  we 
now  proceed  to  tackle  the  title  problem  in  the  following  section. 

Large  Deformation  hi  as t i c - PI  as t ic  Dynamic  B uckling  of  Spherical  Caps 

Prom  Ref.  [10],  we  find  that  the  large  deformation  elastic- 
plastic  response  of  axisymmetric  spherical  caps  under  static 
uniform  loading  is  independent  of  cap  geometric  parameter  X,  a 
factor  used  to  characterize  the  elastic  spherical  cap  behavior, 
and  rather  depends  on  individual  geometric  dimensions.  In  other 
words,  two  different  spherical  caps  with  the  same  X value  have 
the  same  load  carrying  capacity  (p  ) in  elastic  analysis,  but 
may  yield  different  magnitude  of  cl ast i c -pi ast ic  buckling  pres- 
su re  s . 

According  to  this  observation,  two  types  of  spherical  caps 
studied  in  [10],  with  very  much  difference  in  geometrical  di- 
mensions, are  selected  again  for  the  present  analysis.  They 
arc  identified  as  shells  A and  B in  Figs.  7 and  8,  respectively. 

All  geometrical  dimensions  in  shells  A and  B are  fixed  with 
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exception  that  the  thickness  can  be  varied.  By  selecting 
h = 0.26  in.  and  0.0227  in.  respectively  for  shells  A and  B, 
both  caps  yield  the  same  geometrical  parameter  i = 5,  but  have 
different  values  of  the  thickness  to  radius  ratio  (h/R  = 0.0104 
for  shell  A and  0.0048  for  shell  B) . Material  properties  of 
both  caps  given  in  the  figures  are  assumed  to  obey  von  Mises 
initial  yield  condition  with  linear  kinematic  hardening. 

In  Ref.  [6],  two  types  of  dynamic  loading  were  considered 
for  the  buckling  analysis  of  elastic  axisymmetric  spherical 
caps:  uniform  step  loading  with  infinite  duration  and  a right 

triangular  pulse  with  various  time  duration.  Results  obtained 
from  this  reference  reveal  that  pulse  duration  has  a very  signif- 
icant impact  on  the  magnitude  of  the  dynamic  buckling  load,  and 
that  the  step  loading  provides  a more  severe  loading  situation 
than  the  right  triangular  pulse  for  dynamic  analysis.  Based  on 
this  conclusion,  only  uniform  step  loading  with  infinite  duration 
is  considered  for  the  present  analysis. 

As  mentioned  in  SOLUTION  STRATEGY,  a judicious  choice  of 
time  step  At  = 0.1  is  used  for  all  computations  in  this  paper. 

A total  response  time  of  t = 50  is  also  chosen  so  that  if  there 
is  no  sudden  jump  in  the  average  displacement  A or  no  con- 
vergence failure  in  iteration  during  this  period  of  time,  we 
assume  that  no  buckling  occurs  at  this  load.  Convergence 
criteria  are  specified  in  Fig.  5. 

Some  large  deformation  clastic-plastic  dynamic  response 
curves  for  shells  A and  B are  shown  in  Figs.  7 and  8,  respec- 
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tivcly,  which  also  show  shell  A buckles  at  p = 0.26  while 
shell  B at  0.36.  Bor  comparison  purposes,  some  data  of  interest 
are  also  listed  here:  both  caps  have  the  same  elastic  dynamic 
buckling  load  p = 0.46  [6],  while  their  elastic-plastic  static 
buckling  values  are  0.31  for  shell  A and  0.50  for  shell  B [10]. 

According  to  these  results,  some  observations  may  be  in 
order.  First,  buckling  loads  for  spherical  caps  of  elastic 
material  are  related  to  their  geometric  parameter  X.  Secondly, 
plastic  yielding  plays  a significant  role  in  reducing  the 
buckling  pressure  of  spherical  caps;  the  influence  of  plastic 
yielding  increases  with  increase  of  the  thickness  to  radius 
ratio.  For  example,  plastic  yielding  reduces  the  dynamic 
buckling  pressure  pcj.  for  shell  B (h/R  = 0.0048)  from  0.46  to 
0.36,  a 22%  reduction  in  buckling  capacity.  On  the  other  hand, 
the  reduction  increases  to  43‘o  (from  0.46  to  0.26)  for  shell  A 
fh/R  = 0.0104).  This  indicates  that  the  geometric  parameter  A 
becomes  meaningless  in  elastic-plastic  analyses.  Another 
important  observation  is  that  the  dynamic  effect  also  reduces 
the  buckling  pressure  of  spherical  caps.  This  is  evidenced  by 
the  fact  that  dynamic  effect  cuts  down  elastic-plastic  buckling 
loads  for  shell  A from  pc r = 0.31  to  0.26,  for  shell  B from  0.50 
to  0.36. 

From  dynamic  response  curves  in  Figs.  7 and  8,  it  is  found 

that  in  the  load  range  above  the  buckling  pressure  a larger  load 

level  requires  a shorter  response  time  to  reach  A or  to 

max 

generate  convergence  failure  in  iteration.  This  seems  to  be 
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true  for  all  cases  treated  herein  except  for  one  situation 
which  will  be  discussed  when  the  case  is  encountered. 

Because  of  having  a greater  thickness  to  radius  ratio  and 
hence  being  more  sensitive  to  plastic  deformation,  shell  A is 
selected  to  examine  the  effect  of  initial  imperfections  on  the 
cap  elastic-plastic  dynamic  buckling  pressure.  On  the  other 
hand,  shell  B is  studied  again  to  establish  more  evidence  on 
how  the  effect  of  plastic  yielding  on  the  buckling  pressure  is 
influenced  by  the  thickness  to  radius  ratio. 

By  setting  h = 0.0101  in.  and  0.0057  in.  for  shell  B of 
Fig.  8,  two  caps  are  obtained,  one  having  A = 7.5  and  h/R  = 0.0021 
and  the  other  A = 10  and  h/R  = 0.0012.  Large  deformation  elastic- 
plastic  dynamic  response  curves  displayed  in  Figs.  9 and  10 
show  that  buckling  loads  are  read  as  p = 0.43  for  A = 7.5  and 
0.46  for  A = 10,  while  their  corresponding  elastic  dynamic 
values  are  0.44  and  0.49,  respectively  [10].  Reduction  in 
buckling  loads  because  of  plastic  yielding  in  these  two  cases 
may  be  neglected.  Insensitivity  to  plastic  deformation  in  these 
situations  may  be  attributed  to  the  fact  that  as  the  thickness 
to  radius  ratio  becomes  smaller,  the  membrane  effect  tends  to 
be  more  dominant.  At  the  time  when  the  cap  buckles,  the  de- 
velopment of  plastic  deformation  may  be  still  very  limited  or 
of  minor  importance. 

It  is  noted  that  Fig.  9 provides  the  first  and  only  excep- 
tion that  a smaller  load  magnitude  takes  a shorter  response 
time  to  reach  A or  to  yield  convergence  failure  in  itera- 
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tion,  although  the  difference  between  the  maximum  response 
times  for  p = 0.43  and  0.44  is  negligible. 
l Now,  let  us  examine  the  combined  effect  of  nonlinearities 

in  material  and  geometry  and  initial  imperfections  on  the 
dynamic  buckling  of  spherical  caps.  As  already  mentioned 
earlier,  shell  A of  A = 5,  whose  geometric  and  material  pro- 
perties are  given  in  Fig.  7,  is  selected  for  this  analysis. 

Three  values  of  initial  imperfections  BL^/h  = 0.1,  0.5  and  1.0 
imposing  on  the  cap  domain  of  shell  A are  considered.  Dynamic 
response  curves  of  these  three  cases  are  shown  separately  in 
Figs.  11,  12  and  13. 

Buckling  loads  read  from  these  figures  are  pcr  = 0.25, 

0.17  and  0.13  for  BL^/h  = 0.1,  0.5  and  1.0,  respectively.  These 
values  together  with  those  of  static-elastic,  dynamic-elastic 
and  static  elastic-plastic  [6,10]  are  listed  in  Table  1 and 
plotted  in  Fig.  14.  Also  recorded  in  Table  1 for  comparison 
purposes  are  some  buckling  loads  obtained  in  Fig.  8 and  Ref. 

[10]  for  shell  B (A=5). 


Table  1.  Buckling  loads  for  two  clamped  spherical  caps  (shell 
A,  A = 5,  Fig.  7;  shell  B,  A = 5,  Fig.  8) 


w . 

Elastic 

Elastic- 

plastic 

iO 

shells 

A and  B 

shel  1 

A 

shel  1 

15 
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Stat ic 
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Dynamic 
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Static 

[10] 

Dynamic 

Static 
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Dynamic 

0 

. 64 

. 46 

. 31 

. 26 

. 50 

. 36 

. l 

. 54 

. 39 

. 28 

. 25 

. 43 

- 

. 5 

. 32 

.28 

. 20 

.17 

. 27 

. 

1.0  .15 


.12  .13 
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Results  exhibited  in  big.  14  show  that  imperfections  do 
indeed  have  a significant  impact  in  reducing  buckling  capacity 
l for  all  cases  considered  here,  although  they  are  influenced  in 

a different  manner.  The  difference  in  buckling  load  reduction 
may  be  examined  from  two  different  bases.  First,  if  the  ex- 
amination is  focused  on  the  same  material  property  (elastic  or 
clastic-plastic),  the  rate  of  change  of  the  buckling  pressure 
Pcr  with  respect  to  the  imperfection  W^/h  is  smaller  for  the 
dynamic  response  than  for  the  static.  As  a result,  the  reduc- 
tion of  buckling  loads  due  to  dynamic  effect  decreases  with 

* 

increase  of  the  imperfection  magnitude.  As  the  imperfection 
exceeds  a certain  magnitude  ( W . q / h > 0.75  and  0.85  for  the 
clastic  and  elastic-plastic  material,  respectively),  the 
dynamic  buckling  loads  may  be  even  higher  than  the  correspond- 
ing static  buckling  values.  This  finding  suggests  that  the 
statement  --  the  dynamic  effect  introduced  by  a time-dependent 
load  can  reduce  the  load  carrying  capacity  of  a structure  [22] 

--  is  true  for  the  perfect  spherical  cap,  but  not  necessarily 
true  for  the  cap  with  initial  imperfections. 

Secondly,  if  the  examination  is  centered  on  the  same 
loading  situation  (static  or  dynamic),  results  in  Fig.  14 
clearly  demonstrate  that  both  plastic  yielding  and  initial 
imperfection  have  the  same  effect  of  lowering  the  shell  load 
carrying  capacity.  Figure  14  also  shows  that  the  imperfection 
* 

This  reduction  is  referring  to  the  difference  between  static 
and  dynamic  buckling  loads  under  the  same  material  property  and 
the  same  imperfection  magnitude  IV.  ^/h  (see  Fig.  14). 


▲ 
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is  a more  dominant  factor,  since  the  influence  of  plastic 
yielding  diminishes  as  the  imperfection  magnitude  increases. 

The  informations  discussed  here  are  quite  useful  for 
spherical  cap  designs.  A valid  design  of  these  structures 
should  reflect  the  actual  material  property,  consider  both 
static  and  dynamic  responses,  and  also  take  into  account  the 
estimated  magnitude  of  imperfections  in  light  of  manufacturing 
cons iderat ions . 

CONCLUSIONS 

The  objective  of  this  paper  is  to  obtain  large  deformation 
clastic-plastic  dynamic  buckling  loads  of  axisymmetric  spherical 
caps  with  initial  imperfections.  The  problem  formulation  is 
based  on  governing  differential  equations  of  motion,  treating 
the  plastic  deformation  as  effective  plastic  loads  which  are 
combined  together  with  actual  applied  forces. 

Liquations  of  motion  are  converted  into  a discrete  system 
of  equations  by  replacing  all  spacewise  derivatives  with  central 
finite  differences  and  the  second  time  derivative  with  a 
lloubolt's  third-order  backwards  difference  expression  [19]. 
Resulting  nonlinear  algebraic  equations  are  then  solved  by  a 
step-by-step  time  integration  scheme  for  displacements.  At 
each  step  of  integration,  the  solution  procedure  is  divided  into 
two  major  loops,  namely,  elastic  solution  and  material  pro- 
perty loops. 

In  the  elastic  solution  loop,  all  material  properties  are 


held  constunt  and  nonlinear  equations  arc  solved  by  the  non- 
linear relaxation  technique  |12|.  In  the  material  property 
loop,  material  properties  are  updated  to  correspond  the  new 
computed  state  of  stress  to  the  specified  stress - stra in  rela- 
tion. The  procedure  consisting  these  two  loops  is  repeated 
until  governing  equations  and  the  specified  stress  - strain  re- 
lation are  satisfied  at  all  points  over  the  shell  surface  and 
throughout  the  thickness.  The  solution  procedure  discussed 
here  is  sketched  in  Fig.  5 for  references. 

The  plasticity  relation  adopted  herein  is  an  incremental 
flow  theory  which  is  based  on  von  Mises  initial  yield  condition 
and  the  Prager- Z iegler  kinematic  hardening  rule.  It  is  noted 
that  this  theory  takes  into  account  the  Bauschinger  effect  and 
that  the  computer  program  developed  for  the  present  analysis  is 
equipped  to  handle  the  material  of  elast ic-perfectly  plastic, 
linear  and  nonlinear  hardening  behavior.  Nine  thickness  points 
and  a Simpson  rule  are  used  for  all  numerical  calculations  of 
plastic  deformation  in  this  paper.  A dimple  type  of  imperfec- 
tions is  selected  whicli  provides  a quite  aaequite  description 
of  the  local  nature  of  spherical  shells. 

To  verify  the  validity  of  the  theory  and  the  solution 
method  employed,  the  procedure  is  applied  to  obtain  the  large 
deformation  clastic-plastic  dynamic  response  of  a spherical 
cap  (Fig.  6).  The  present  solution  along  with  those  obtained 
by  Nagarajan  and  Popov  [20]  and  Bathe,  Ramm  and  Wilson  [21] 
arc  displayed  in  Fig.  6.  Their  solutions  are  based  on  a finite 


element  formulation,  isotropic  strain  hardening  and  a Newmark 
time  integration  scheme.  Despite  differences  in  the  theory 
and  solution  methods  utilized,  the  comparison  among  these  three 
sets  of  results  is  remarkably  good . 

The  solution  procedure  is  then  applied  to  solve  the  title 
problem;  two  spherical  cap  models  shown  in  Figs.  7 and  8 are 
chosen  for  the  present  analysis.  The  loading  exerted  on  the 
cap  surface  is  assumed  to  be  a uniform  step  pressure  of  infinite 
duration.  The  results,  which  are  given  and  discussed  in  great 
details  in  the  previous  section,  may  be  summarized  as  follows: 

(1)  For  axisymmetric  spherical  caps  of  elastic  material, 
their  static  and  dynamic  buckling  loads  (prr)  are  function  of 
cap  geometric  parameter  X.  However,  for  caps  of  elastic-plas- 
tic material,  these  values  are  independent  of  X.  This  means 
that  two  spherical  caps,  which  have  the  same  X value  but  dif- 
ferent geometric  dimensions,  may  yield  the  same  magnitude  of 
buckling  load  (Pcr)  if  their  materials  are  elastic;  however, 
their  pcr  values  may  be  different  if  their  materials  are  elas- 
tic-plastic. For  example,  Table  1 shows  that  if  both  shells  A 
and  B (both  have  X = 5)  are  made  of  the  elastic  material,  their 
static  buckling  loads  are  of  the  same  magnitude  (Pcr  = 0.46); 
both  also  have  the  same  dynamic  buckling  value  (Pcr  = 0.46). 

But  when  their  materials  are  assumed  to  be  elastic-plastic, 

the  static  buckling  loads  are  0.31  and  0.50,  and  the  dynamic 
values  are  0.26  and  0.36  for  shells  A and  B,  respectively. 

(2)  Plastic  yielding  plays  a significant  role  in  reducing 
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the  buckling  pressure  of  spherical  caps.  The  influence  of 
plastic  yielding  increases  with  increase  of  the  thickness  to 
radius  ratio.  This  observation  may  be  illustrated  from  Table 
1 that  plastic  yielding  reduces  the  dynamic  buckling  pressure 
Pcr  for  shell  B (h/r  = 0.0048)  from  0.46  to  0.36,  a 22%  re- 
duction in  buckling  capacity.  On  the  other  hand,  this  reduc- 
tion increases  to  43%  (from  0.46  to  0.26)  for  shell  A (h/R  = 0.0104). 

(3)  Tor  perfect  spherical  caps,  dynamic  effect  has  an 
influence  of  lowering  their  load  carrying  capacity.  This  can 

be  seen  from  Table  1 that  dynamic  effect  cuts  down  their  elastic- 
plastic  buckling  loads  (p  ) from  0.31  to  0.26  and  from  0.50 
to  0.36  for  shells  A and  B,  respectively. 

(4)  For  imperfect  spherical  caps,  dynamic  effect  on  their 
load  carrying  capacity  depends  on  the  magnitude  of  imperfections. 

Tor  example,  the  elastic  rc'sults  of  shell  A ( A = 5 ) displayed  in 
Tig.  14  demonstrate  that  when  imperfection  W^/h  is  less  than 
0.75,  its  dynamic  pcr  values  are  smaller  than  its  static  values. 

On  the  other  hand,  when  W.^/h  exceeds  0.75,  the  situation  is 
completely  reversed.  A similar  situation  is  also  observed  for 
shel 1 A of  elastic-plastic  material. 

Finally,  it  may  be  concluded  that  a realistic  design  of 
spherical  caps  should  be  based  on  the  actual  material  property, 
consider  both  static  and  dynamic  responses,  and  take  into 
account  the  estimated  magnitude  of  imperfections  in  light  of 
ma nu factur ing  considerations. 
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Fig.  5 — Iteration  procedure  for  large  deformation  elastic - 
plastic  dynamic  problems. 


Fig.  6.  Large  displacement  dynamic  elastic -plastic  analysis 
of  spherical  cap  under  uniform  external  pressure  of 
600  psi. 
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Fig.  7 Dynamic  response  for  a clamped,  elastic-plastic  spherical  cap  (A  = 5)  under  uniform  step  pressure 
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Fig.  9.  Dynamic  response  for  a clamped,  elastic-plastic  spherical  cap  ( X = 7. 5 ) under  uniform  step  pressure. 
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Fig.  11.  Dynamic  response  fora  clamped,  elastic- plastic  spherical  cap  (X=5)with  initial  imperfection  Wj/h=  0.1  and 
under  uniform  step  pressure  (Shell  A). 


Fig.  12.  Dynamic  response  for  a clamped,  elastic- plastic  spherical  cap  (X=5)  with  initial  imperfection  Wj/h  = 
and  under  uniform  step  prpssure.  (Shell  A). 


Non-dimensional  Time,  T 

Fig.  13.  Dynamic  response  for  a clamped,  elastic- plastic  spherical  cap  (A,=  5)  with  initial  imperfection  Wio=1.0 
and  under  uniform  step  pressure.  (Shell  A). 
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Fig.  14.  Buckling  load  vs.  initial  imperfection  for  an  axisymmetric  spherical  cap  (Shell  A,A=5,  Fig.  7). 
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